Roughening and preroughening transitions in crystal surfaces with double-height steps 
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We investigate phase transitions in a solid-on-solid model where double-height steps as well as 
single-height steps are allowed. Without the double-height steps, repulsive interactions between 
up-up or down-down step pairs give rise to a disordered flat phase. When the double-height steps 
are allowed, two single-height steps can merge into a double- height step (step doubling). We find 
that the step doubling reduces repulsive interaction strength between single-height steps and that 
the disordered flat phase is suppressed. As a control parameter a step doubling energy is introduced, 
which is assigned to each step doubling vertex. From transfer matrix type finite-size-scaling studies 
of interface free energies, we obtain the phase diagram in the parameter space of the step energy, 
the interaction energy, and the step doubling energy. 



Much attention has been paid to the phase transitions 
in crystal surfaces since they show rich critical phenom- 
ena. The interplay between roughening and reconstruc- 
tion results in interesting phases, such as a disordered 
flat (DOF) phase, as well as flat and rough phases [Q. 
In the DOF phase the surface is filled with macroscopic 
amount of steps which are disordered positionally but 
have up-down order. Several solid-on-solid (SOS) type 
models have been studied, which reveals that the DOF 
phase is stabilized by the repulsive step-step interac- 
tions l^-Q or by specific topological properties of sur- 
faces, e.g., Si(OOl) g. 

The SOS type model studies have been done in cases 
where the nearest-neighbor (NN) height difference, A/i, is 
restricted to be equal to or less than 1 in units of the lat- 
tice constant. However, in real crystals there also appear 
steps with |A/i| > 1. For example, double-height steps 
on W(430) become more favorable than single-height 
steps at high temperatures since they have lower kink 
energy g. In this paper we investigate the phase tran- 
sitions in crystal surfaces in the presence of the double- 
height steps with I A/i| = 2, especially focusing on the sta- 
bility of the DOF phase. We study a generalized version 
of the restricted solid-on-solid (RSOS) model on a square 
lattice with the Hamiltonian given in Eq. (||). We study 
the model under the periodic and anti-periodic boundary 
conditions, from which various interface free energies are 
defined. The interface free energy is calculated from nu- 
merical diagonalizations of the transfer matrix, and the 
phase diagram is obtained by analyzing their finite-size- 
scaling (FSS) properties. 

In the RSOS model the surface is described by integer- 
valued heights hr at each site r = (n, m) on a square 
lattice. (The lattice constant in the z direction is set to 
1.) Only the single- height step (S step) with |A/i| = 1 
is allowed. It was found that the RSOS model with NN 
and next-nearest-neighbor (NNN) interactions between 
height displays the DOF phase when the NNN coupling 
strength is large enough . The NNN coupling accounts 



for the repulsive interactions between parallel (up-up or 
down-down) step pairs. Parallel step pairs cost more en- 
ergy than anti-parallel (up-down or down-up) step pairs. 

The double-height step (D step) is incorporated into 
the RSOS model by relaxing the restriction on the NN 
height difference to |A/i| = 0,1,2. We only consider 
quadratic NN and NNN interactions between heights 
since they are sufficient to describe the key feature of 
the phase transitions. The total Hamiltonian is written 
as 



(1) 



(r,r') 



(r.r") 



where () and () denote the pair of NN and NNN sites. 
With this Hamiltonian, a D step costs more energy than 
two separate S steps by an amount of 2K + AK per unit 
length. Even though the D steps are energetically unfa- 
vorable, we will show that their effect is not negligible. 
We also consider a step-doubling energy Ed to study the 
effect of the step doubling. It is assigned to each vertex 
where two S steps merge into a D step (see Fig. |l|). The 
electronic state at step edges may be different from that 
at a flat surface, which contributes to the step energy. 
When two S steps merge into a D step, the electronic 
state near the vertex may be changed. The change leads 
to an additional energy cost, which is reflected by Eo. 
When Ed is positive (negative), it suppresses (enhances) 
the step doubling. The Hamiltonian including Hq and 
the step-doubling energy is then given by 



H = Ho + EdNl 



(2) 



where Njj is the total number of step-doubling ver- 
tices. (For a notational convenience the energy is mea- 
sured in unit of fc^T.) The model with the Hamiltonian 
Eq. (||) with Ed — and with the restriction \Ah\ = 0, 1 
will be referred to as the RS0S3 model, and the model 
with the Hamiltonian Eq. (|^) and with \Ah\ = 0, 1, 2 will 
be referred to as the RS0S5 model. 
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FIG. 1. Step doubling vertices. 

In a continuum description phase transitions in crystal 
surfaces are described by the sine-Gordon model 



H 



(3) 



where (j){r) G (— oo,cx3) is a real- valued local average 
height field, Kq the stiffness constant, and Uq the fugac- 
ity of q-charge 0. In the renormalization group sense ui 
is irrelevant at high temperatures where the model renor- 
malizes to the Gaussian model with a renormahzed stiff- 
ness constant Kq < f describing the rough phase. As 
temperature decreases, ui becomes relevant at a rough- 
ening transition temperature. There appear two kinds 
of low temperature phases depending on the sign of ui: 
For positive ui the Hamiltonian favors an integer aver- 
age height and hence the surface is flat. For a negative 
Ml it favors a half-integer average height. Since the mi- 
croscopic height is integer-valued, the surface can take 
the half-integer average height by forming steps with up- 
down order, i.e., the surface is in the DOF phase. As 
temperature decreases further, the sign of ui changes and 
the surface falls into the flat phase. At the roughening 
transition between the rough phase and the flat or DOF 
phase, the renormalized stiffness constant takes the uni- 
versal value of ^. The flat and DOF phases are sepa- 
rated by the preroughening transition characterized by 

Ml = 0. 

The phase boundaries can be obtained using FSS prop- 
erties of the interface free energies. Consider the model 
on a finite N x M square lattice rotated by 45° under the 
various boundary conditions (BG's): The periodic BC, 
h{n -j- iV, m) = h{n, m) -\- a with integer a, and the anti- 
periodic BC, h{n + N, m) = —h{n, m) + a (mod 2) with 
a — and 1. They will be referred to as (±, a) BC's (the 
upper (lower) sign for the (anti-)periodic BC's). The 
free energy is obtained from the largest eigenvalue of the 
transfer matrix. Detailed description of the transfer ma- 
trix set-up can be found in Ref. (^J^. The boundary 
conditions except for the (+,0) BC induce a frustration 
in the surface. The interface free energy 7/^ is defined as 
the excess free energy per unit length under the k BC 
with K = (±, a) from that under the (+, 0) BC: 



M Z, 



( + .0) 



(4) 



with Zf^ the partition function satisfying the k-BC. 

The interface free energies have characteristic FSS 
properties in each phase. In the rough phase they 



show the universal I /N scaling in the semi- infinite limit 

M ^ 00 as 



2 N 



(5) 
(6) 



where Kq < f is the renormalized stiffness constant of 
the Gaussian model and C, is the aspect ratio of the lattice 
constants in the horizontal and vertical directions 
In the flat phase ??(+,a) s^nd ili-.i) ^re finite because at 
least one step is induced under the (-I-, a) and (— , 1) BC's, 
while ??(-,o) is exponentially small in N since the (— ,0) 
BC may not induce any steps |||. In the DOF phase 
the (— , 1) BC does not induce any frustration in the step 
up-down order, but the (+,a) and (—,0) BC's do. So 
is exponentially small in N, and ^(+^a) ^-nd '7(-,o) 
are finite From these FSS properties the roughening 
points can be estimated from 



4iV' 



(7) 



where the universal value of Kg = ^ at the roughening 
transition is used in Eq. (|^). The preroughening points 
between the flat and the DOF phase can be estimated 
from the crossing behaviors of iV77(_ 0) or -^^(-,1): which 
converges to zero in one phase and diverges to infinity in 
the other phase as N grows. 

The estimation of transition points using the interface 
free energies suffers from slow convergence due to correc- 
tions to the scaling. They may smooth out the crossing 
behaviors of iVry^. 0) and A^77(__i) at the preroughening 
transitions for small N . But one can safely cancel out 
leading corrections to scaling by taking the ratio or the 
difference of them, which can be seen as follows. Consider 
the lattice version of the continuum model in Eq. (||) . It 
is obvious, using the transformation (j) (j) — 1/2, that 
the model under the (—,0) BC is the same as that un- 
der the (— , 1) BC with Uq replaced by —Uq for odd q. It 
yields the relation 

?7(_,0)(W1,M2,W3, ■ • •) = »7(-4)(-'"i'"2,-U3, ■ • •) ■ (8) 

So if one neglects all higher order contributions from 
Mq>3, the location of Mi = is found from the condi- 
tion ??(-,o) ~ ^(-,1) = or i? = 1 with 



(9) 



It is not influenced by correction-to-scalings from U2. 
Therefore the relation R ~ 1 can be used to get the 
Ml = point more accurately. One can easily see that 
i? > 1 for negative mi and R < 1 for positive mi. It 
approaches 1 in the rough phase and at the prerough- 
ening transition points, diverges in the DOF phase, and 
vanishes in the flat phase as ^ 00. 
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FIG. 2. The ratio i? defined in the text for the RS0S3 
model (a) and the RS0S5 model (b) along the line L = ^K. 
The rough, DOF, and flat phases are denoted by R, DOF, 
and F, respectively. (□ for iV = 4, o for = 6, A for = 8, 
and V for TV = 10.) 

In the RS0S3 model the exact point with mi = is 
known along the line i = 0. It is called the self-dual 
point and is located a,i K = Ksd — ln[|(\/5 + 1)]. From 
numerical studies of the RS0S3 model transfer matrix, 
we could obtain the exact value of Ksd with error less 
than 10~^^ by solving R = 1 even with small system size 
TV = 4, which indicates that i? is a useful quantity to de- 
termine the preroughening transition points accurately. 
It will be used in the analysis of the RS0S5 model. 

We first consider the RS0S5 model in a special case 
of Ed = and compare its phase diagram with that of 
the RS0S3 model to have insight into the role of the D 
step. At low temperatures the D step is unfavorable due 
to larger free energy cost than the S step. So the nature 
of the low temperature phase in the RS0S5 model is not 
different from that in the RS0S3 model, i.e., the flat 
phase. At high temperatures, the surface is in the rough 
phase in the RS0S3 model. Since the rough phase is crit- 
ical and there is no characteristic length scale, there will 
be no difference between S and D steps. So the RS0S5 
model will also have the rough phase as a high tempera- 
ture phase. There is significant difference in intermediate 
temperature range, where the repulsive step interactions 
stabilize the DOF phase in the RS0S3 model. With- 
out the D steps the parallel steps have less meandering 
entropy than anti-parallel ones. It is energetically unfa- 
vorable for parallel steps to approach each other closer 
than the interaction range while anti-parallel steps can 
approach each other at will However, if one allows 
the D step, two parallel S steps can approach each other 
and form a D step without the interaction energy cost. 
Provided that the energy cost of the D step is not too 
high, the presence of the D step reduces repulsive inter- 



action strength effectively and enhances the meandering 
entropy of parallel steps. Then it will suppress the DOF 
phase. 

To see such effects quantitatively, we calculate the ra- 
tio R for the RS0S3 model and the RS0S5 model with 
Ed ~ along a line L = 5K (see Fig. ||). The strip 
width for the transfer matrix is TV = 4, 6, 8, and 10 for the 
RS0S3 model and TV = 4, 6, and 8 for the RS0S5 model. 
The RS0S3 model displays the roughening and the pre- 
roughening transitions along the line L = 5K, which is 
manifest in Fig. ^(a). There are three regions where A'' 
dependence of R is distinct with each other. The sur- 
face is in the rough phase with negative ui in the small 
L (high temperature) region, where R approaches 1 from 
above. And the surface is in the DOF (flat) phase for the 
intermediate (large) L region, where R grows (vanishes). 
The roughening and preroughening transition points are 
estimated from Eq. (^ and R = 1 with R in Eq. (||), re- 
spectively, which are represented by broken vertical lines. 

The situation changes qualitatively in the RS0S5 
model. As can be seen in Fig. ||(b), R is always less 
than 1, and there are only two regions with distinct N 
dependence of R. In the small L region R approach 1 
from below, and in the large L region R vanishes as N 
increases. They correspond to the rough phase with pos- 
itive ui and the flat phase, respectively. The roughening 
transition point is estimated from Eq. (R) and represented 
by the broken vertical line. It shows that the DOF phase 
is suppressed in the presence of the D step. We have also 
checked that R is always less than 1 {ui > 0) and the 
DOF phase does not appear at any values of K and L in 
the RS0S5 model with Ed = 0. 
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(a) (b) 
FIG. 3. Comparison of configurations of two parallel steps 

which cost the interaction energy (a) and from a D step (b). 

S (D) steps are denoted by single (double) arrows. 

We can argue the reason why the DOF phase disap- 
pear in the presence of the D step as follows. Consider 
two parallel S steps merging at a vertex. If the D step is 
not allowed, the possible vertex configuration is as shown 
in Fig. ||(a) and the energy cost for such configuration is 
2K + AL. On the other hand, if the D steps is allowed, 
the step doubling may occur in two ways as shown in 
Fig. 1(b) with the energy cost 3K + 5L. Though the 
step doubling costs more energy [K + L) , entropic con- 
tribution of the step doubling (— In 2) may lower the free 
energy of parallel steps below than the value without the 
step doubling. Our numerical results above show that the 
step doubling suppresses the DOF phase entirely in the 
Ed = case. In our model a D step costs more energy 
than two separate S steps. The two energy scales may 
be comparable to each other in a more realistic model. 
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FIG. 4. The ratio R for the RS0S5 model with the step FIG. 5. The phase diagram of the RS0S5 model in the 

doubling energy e'^" = 0.2 along the line L = 5K. (□ for {K, L, e~'^'=') parameter space. The roughening transition 
N = i, o for A'^ = 6, and A for TV — 8.) hnes are denoted by broken lines and the preroughening tran- 

sition lines by solid lines. 



From the above arguments, one finds that the step dou- 
bUng plays an important role in phase transitions. So we 
introduce a new term EdNd in Eq. (|^) with the step- 
doubling energy Eb and study the phase diagram in the 
parameter space {K,L,Ed)- When Ed < 0.0 (> 0.0), 
the step doubling is favored (suppressed). One can easily 
expect that the DOF phase does not appear for negative 
Ed- 

For positive Ed the step doubling is suppressed and the 
effect of the step interaction becomes important. So we 
expect there appears the DOF phase in the positive Ed 
side of the parameter space. In Fig. ^ we show the ratio 
i? for e"-^° =0.2 and along the line L = 5K. Though 
the convergence is not good, compared with Fig. ^(a), 
one can identify three regions as the rough, DOF, and 
flat phases from the N dependence of R. The roughen- 
ing point between the rough phase and the DOF phase 
is estimated using Eq. and the preroughening point 
using R — 1 for = 8. They are denoted by broken 
vertical lines. 

We obtain the phase diagram in the whole parameter 
space using the conditions = ^'^^ roughen- 

ing transition boundary and R = 1 for the preroughen- 
ing transition boundary. It is obtained for strip width 
iV = 4,6, and 8. Since the maximum TV we can handle 
is small, the convergence of the phase boundary is poor 
especially as one approaches e~'®° = 0. But there is no 
qualitative change in shape. So we only present the phase 
diagram obtained from TV = 8 in Fig. a. The region un- 
der the surface represented by broken lines corresponds 
to the rough phase. The DOF phase is bounded by the 
surfaces of broken lines and solid lines. The region above 
the surfaces corresponds to the flat phase. One should 
notice that there is a critical value of Ed, approximately 
0.071, smaller than which the DOF phase does not ap- 
pear. 



In summary, we have studied the phase transitions in 
the RS0S5 model with the Hamiltonian in Eq. (|^) with 
D steps as well as S steps. We have found that the D step, 
which has not been considered in previous works, plays 
an important role in phase transitions in crystal surfaces. 
The presence of the D step reduces the strength of the 
repulsive interaction between parallel steps through the 
step doubling, and hence suppresses the DOF phase. We 
also found that the step-doubling energy is an important 
quantity which characterizes a surface upon the rough- 
ening. 
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